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ABSTRACT 

The methanol molecule CH3OH has a complex microwave spectrum with a large number of 
very strong lines. This spectrum includes purely rotational transitions as well as transitions with 
contributions of the internal degree of freedom associated with the hindered rotation of the OH 
group. The latter takes place due to the tunneling of hydrogen through the potential barriers 
between three equivalent potential minima. Such transitions are highly sensitive to changes in 
the electron-to-proton mass ratio, /1 = m c /m p , and have different responses to /i-variations. 
The highest sensitivity is found for the mixed rotation-tunneling transitions at low frequencies. 
Observing methanol lines provides more stringent limits on the hypothetical variation of [i than 
ammonia observation with the same velocity resolution. We show that the best quality radio 
astronomical data on methanol maser lines constrain the variability of /x in the Milky Way at the 
level of \A[i/[i\ < 28 x 1CP 9 (lcr) which is in line with the previously obtained ammonia result, 
\Afi/fi\ < 29 x 10" 9 (lcr). This estimate can be further improved if the rest frequencies of the 
CH3OH microwave lines will be measured more accurately. 



Subject headings: molecular data 
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Introduction 



techniques: radial velocities — ISM: molecules — dark energy 



The hypothetical variation of the dimensionless 
physical constant p, - the electron-to-proton mass 
ratio - can be probed through the spectral ob- 
servations of certain molecular transitions which 
are particularly sensitive to changes in [i. The 
corresponding sensitivity coefficients Q M of differ- 
ent molecular transitions relevant to astrophysics 
were calculated at first for H2 (Varshalovich & 
Levshakov 1993), and later on for OH (Darling 



2003; Chengalur & Kanekar 2003), 15 ND 3 (van 
Veldhoven et al. 2004), H 2 (Ubachs et al. 2007), 
NH 3 (Flambaum & Kozlov 2007), OH and CH 
(Kozlov 2009), NH 2 D and ND 2 H (Kozlov et al. 
2010), H 3 0+ (Kozlov & Levshakov 2011), CH 3 OH 
(Jansen et al. 2011), and H 3 0+, H 2 DO+, HD 2 0+, 
and D 3 + (Kozlov et al. 2011). Among them am- 
monia, NH 3 , is actively used in extragalactic and 
galactic observations of dense molecular clouds. In 
extragalactic observations the most stringent lim- 
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its on the mass ratio, A/z//x = (//space - Wab)/Wab, 
were obtained by Henkel et al. (2009) at redshift 
z = 0.89 and by Kanekar (2011) at z = 0.69: 
\Afi/fj,\ < 1400 ppb (3a-), and |A/VmI < 400 ppb 
(3<t), respectively (lppb = 10 -9 ). Observations 
in the Milky Way have shown, however, a ten- 
tative signal An/n — 26 ± 3 ppb (Levshakov et 
al. 2010a, b), but its nature remains unclear. To 
distinguish whether this is a real signal or an arte- 
fact due to unaccounted systematic effects, addi- 
tional independent measurements involving other 
molecules are required. 

One of such molecules — CH3OH — was re- 
cently suggested by Jansen et al. (2011, here- 
after Jll). CH3OH is a widespread interstel- 
lar molecule observed in the Milky Way, external 
galaxies (Herbst & van Dishoeck 2009; Sjouwer- 
man et al. 2010; Martin et al. 2006), and even in 
comets (Bockelee-Morvan et al. 1991). The pur- 
pose of the present paper is to probe the variability 
of /i in the Milky Way using narrow emission lines 
of the methanol masers. 

The microwave spectrum of CH3OH is very rich 
because of the internal rotation of the OH group. 
The basic theory of the non-rigid tops with inter- 
nal rotation was established in the fifties (Lin & 
Swalen 1959; Herschbach 1959) and the main fea- 
tures of the methanol spectrum were explained. 
Later on the theory was refined many times and 
currently there is a very impressive agreement be- 
tween the theory and experiment (Anderson et al. 
1990; Miiller et al. 2006; Xu et al. 2008; Kleiner 
2010). 

2. Effective Hamiltonian and sensitivity 
coefficients 

In this section the sensitivity coefficients Q M of 
the methanol lines are calculated. Our computa- 
tional procedure differs from that described in Jll. 
This allows us to check the values of and their 
validity. 

To calculate the microwave spectrum of CH3OH 
we use an approach from Rabli & Flower (2010, 
hereafter RF), where a simple and convenient form 
of the effective Hamiltonian with six spectroscopic 
constants is suggested. This Hamiltonian is physi- 
cally transparent and sufficiently accurate for cal- 
culations of sensitivity coefficients. All six param- 
eters of the RF model have clear physical meaning 



and their dependence on fi is easily understood 
within the Born-Oppenheimer approximation. 

To show how rotational parameters scale with 
\x, we consider an example of a diatomic molecule 
in its ground vibrational state: 
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where B e = 1/(MRq) corresponds to the equilib- 
rium internuclear distance Rq and a e is the vibra- 
tional correction. It is clear that B e ~ /i, but a e 
has an additional dependence on fi via the vibra- 
tional wave function and, hence, scales as /z 3 / 2 . 
Thus, there is no one-to-one correspondence be- 
tween terms of the effective Hamiltonian and the 
terms of the Born-Oppcnhcimcr perturbation the- 
ory. As a result, B scales as SB /B = Q^Sfj,/ V, 
with Q r w 1. Typically, a e /Bo is of the order 
of 10 and Q r ^ varies from 0.995 for the NO 
molecule to 0.981 for the H-bearing HF molecule. 

We conclude that, in general, the rotational pa- 
rameters A, B, and C scale linearly with fi within 
the uncertainty interval of 1-2%. These vibra- 
tional corrections are of the same order of magni- 
tude as centrifugal corrections considered in Jll. 
Both types of corrections are included below in the 
estimate of the error of the calculated sensitivity 
coefficients. 

The rotational part of the Hamiltonian H TO t 
corresponds to the slightly asymmetric top and in- 
cludes rotational constants A, B, and C {B k, C). 
Here we use the standard convention A > B > C, 
while in RF C > B > A. The hindered rotation is 
described by the Hamiltonian 



H F d + Va h 



■ cos 3a; ) , (2) 



where the kinetic coefficient F is proportional to /j, 
and the electronic potential V3 is independent on 
fi. This model does not include centrifugal distor- 
tions. Interaction of internal rotation with over- 
all rotation is described by a single parameter D, 
which scales linearly with /j, (Lin & Swalen 1959). 
For the rotational degrees of freedom we use the 
basis set of the prolate symmetric top and plane 
waves exp(z77iu;) for the internal rotation. The 
Hamiltonian ([2J mixes waves with m! — m = 3n. 
All relevant matrix elements are tabulated in RF. 
The fitted values of the parameters are also given 
there. 
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After the effective Hamiltonian is formed it is 
diagonalized numerically. Due to the C3 symme- 
try of the Hamiltonian ([2]) the final eigenstates can 
be classified as A-type states and twofold degen- 
erate E-type states. A-type states have definite 
parity p: for A + states p = (— f) J and for A~ 
states p = (— 1) J+1 . Because methanol is close 
to the symmetric top its states are classified with 
an approximate quantum number K, which cor- 
responds to the projection of the angular momen- 
tum J on the axis of the CH3 group. In this study 
we are interested only in the lowest states of the 
internal rotation of A and E symmetry: excita- 
tion of higher states requires kinetic temperatures 
T > 300K, while warm molecular clouds have typ- 
ical temperatures T < I00K. 

The CH3OH transitions with AK — can 
be considered approximately as rotational, where 
the state of the internal motion does not change. 
These transitions have 'normal' sensitivities to \±- 
variation with « 1. Because our model does 
not include centrifugal corrections, must be 
exactly equal to 1. We estimate that these cor- 
rections do not exceed 3% for J < 10. The Jll 
model accounts for these corrections and provides 
the sensitivity coefficients for the AK = tran- 
sitions between 1.00 and 1.03 in agreement with 
our estimates. We note that this deviation from 
unity is of the same order of magnitude as the 
vibrational correction which is not accounted for 
in both the RF and Jll effective Hamiltonians. 
The transitions with AK = ±1, on the contrary, 
lead to the change in the internal motion and have 
sensitivities to /i-variation which vary in a wide 
range. Such transitions are of primary interest for 
our purpose. 

In order to determine the sensitivity coefficients 
Q M for the microwave transitions we first find the 
dependence of the eigenvalues Ei on A/1//1: 



AEi = qi -Z 
A* 



(3) 



where the coefficient individual for each level, 
shows a response of the level Ei to a small change 
of [i {\A[ij [i\ <C 1). This is done by diagonal- 
izing the effective Hamiltonian for three sets of 
the parameters, which correspond to /1 = /io, and 
fj, = [Xq(1 ± e), where e equals to 0.001 or 0.0001 
(in both cases we obtain the same g-factors) . The 
parameters A, B, C, D, and F, all scaling linearly 



with /j,, should then take the values A , A (l ± e), 
etc., where the first value Aq corresponds to the 
fit of the experimental spectrum. The dimension- 
less sensitivity coefficient for the transition 
uj = E up — E\ ow is found through calculation of 
the corresponding q-factor 



where q = q up - q low , and 



Qm = - • 

U! 



(4) 



(5) 



The results of calculations for transitions with 
AK = ±1 are listed in Table [TJ where we also give 
experimental and calculated transition frequen- 
cies. One can see that the accuracy of our model is 
about 1 GHz for all considered transitions. Con- 
sequently, the relative error for the high frequency 
transitions is only a fraction of a percent, but for 
the lowest frequency 6.7 GHz it is 15%. In order 
to improve the accuracy of the Q^-values we use 
experimental frequencies in Eq. ([5]). To check the 
computational procedure we also determined 
for a few transitions with AK = (Qo - loA + 48.4 
GHz, 1 - 2 A+ 96.9 GHz, l x - 2 X A+ 96.1 GHz, 
li - 2iA~ 97.7 GHz, and 4 2 - 5 2 A- 242.2 GHz) 
and for all of them found = 1.00. 

Table [1] also lists the errors of the sensitivity 
coefficients in their last digits which are given in 
parenthesis. These errors were estimated in the 
following way. The effective Hamiltonian from 
RF includes the rotational and tunneling parts as 
well as the interaction between them. The transi- 
tion frequencies and the sensitivity coefficients are 
found from the numerical diagonalization of this 
Hamiltonian followed by the numerical differenti- 
ation of the results in respect to /1. 

Only the low frequency transitions can have 
an enhanced sensitivity to the fi- variation because 
they are of a mixed character with the rotational 
and tunneling contributions to the transition en- 
ergy: 

u> = u r — L} t . (6) 

These two contributions have the following depen- 
dences on fi: 



A/i Aujt 



UJ, 



A/i 



(7) 
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Then, the resultant sensitivity of the transition u 



is given by: 
Auj 



U! J [I 



(8) 



Let us suppose that the sensitivities Q r ^ and 
Q* are known with a relative error e. Then, the 
absolute error of <5 M is equal to: 



AQ, 



£ X 



Q r u uj r + Q*u t 



(9) 



In order to use this expression to estimate the 
errors AQ^ we need to know the decomposition 
© and the sensitivity coefficients <3„ and Q* . As 
we discussed earlier, Q r ^ w 1. We can estimate the 
tunneling sensitivity to be « 2.6 from the semi- 
classical (WKB) approximation (e.g., Kozlov et al. 
2010). The tunneling part of the transition energy 
u>t can be estimated from the model of Hecht & 
Dcnnison (1957). The rotational energy uj r is then 
given by the experimental frequency w and Eq. ([6]) . 

We conservatively estimate the relative error of 
our calculations of Q r ^ and Q 1 ^ to be of about 3%, 
i.e., e — 0.03. This error is associated with the 
missing centrifugal distortion and vibration cor- 
rection which changes the rotational and tunnel- 
ing frequencies by a few percent at the most. The 
resulting errors AQ^ are estimated from Eq. ((9]). 

Table [1] shows that the sensitivity coefficients 
of the mixed rotation-tunneling transitions at low 
frequencies lie in the range between —17(1) and 
+43(3). For comparison, the H2 Lyman and 
Werner transitions have ~ 0.03 (Varshalovich 
& Levshakov 1993), and the ammonia inversion 
transition has Q M = 4.46 (Flambaum & Kozlov 
2007). This means that methanol is almost 1000 
times more sensitive to the change in fi than 
molecular hydrogen and about 10 times more sen- 
sitive than ammonia. Besides, CH3OH has many 
lines with different sensitivities of both signs which 
allows us to estimate A/i//z from observations of 
only this one molecule. This is the advantage over, 
e.g., ammonia NH3 which has the same sensitiv- 
ity coefficients for the inversion transitions (1,1), 
(2,2), etc., and, hence, requires some rotational 
transitions of other molecule as a reference in or- 
der to trace A/if/i. 

Now we can compare our values with the 
sensitivity coefficients derived in Jll. We note 



that jj, was defined in Jll as m p /m e and, hence, 
one expects = —Ku. This comparison is shown 
graphically in Fig. [TJ The error bars in this fig- 
ure mark the la uncertainties. In total we have 
35 common computations of the sensitivity coef- 
ficients. The theoretical calculations of their val- 
ues show a good concordance except for the six 
points at 9.936, 37.703, 38.293, 38.453, 86.615, and 
86.902 GHz where the discrepancy + ex- 
ceeds the 2a level (in Fig.[U the corresponding 
values are -14, 5.1, 12.1, 12.1, 5.9, and 5.9). The 
reason for such deviations is not clear. The highest 
offsets of 2.8a and 2.7cr are found for the J — 5 — 6 
transitions at 38.293 and 38.452 GHz, respectively, 
where the RF model should have a sufficiently ac- 
curate result. On the other hand, the errors of the 
values were taken in Jll to be 5% if \K^\ > 1 
(relative error) or 0.05 if \K^\ < 1 (absolute error). 
However, analysis based on Eq. (j9]) shows that the 
relative errors noticeably differ from line to line. 

3. Observational constraints on A/i//x 

The agreement between the values determined 
with independent methods gives us confidence 
that, given the high values with different signs, 
methanol is excellent for testing the electron-to- 
proton mass ratio. Below we consider such a test 
based on high angular and high spectral resolution 
observations of CH3OH maser lines. 

Equations (jlj and ([3]) show that for a given 
transition from Table [1] LOi, with the sensitivity co- 
efficient Qi, the expected frequency shift, Auji/uji, 
due to a change in fi is given by 



Ac 



A/Li 

A* 



(10) 



Then the value of A/1//1 can be estimated from 
two transitions with different sensitivity coeffi- 
cients Qi and Qy. 



A/i 



Vj ~ Vi 
c(Qi -Qj)' 



(11) 



where Vj and Vi are the apparent radial velocities 
of the corresponding CH3OH transitions and c is 
the speed of light. 

Interstellar CH3OH lines were widely observed 
in the last two decades and we can obtain some 
preliminary estimates of A/i//i using the published 
data. In particular, the maser CH3OH emission is 
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Table 1: Numerical calculation of the Q- factors for the low frequency mixed rotation-tunneling transitions 
(AK = ±1) in methanol. The rest frequencies are taken from Miiller et al. (2004)''' except for those marked 
by the asterisk which are from Lovas (2004). The la uncertainties of the least significant figure of the rest 
frequencies and Q M -values are given in parenthesis. 
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Note. — tThe uncertainties of the first three transitions at 6.668, 9.936, and 12.178 GHz given in Miiller et al. (2004) 
correspond to the 2a errors in accord with Breckcnridgc & Kukolich (1995). 



of a special interest here since maser lines are nar- 
row and the error of the line center measurement 
is correspondingly low. 

The CH 3 OH molecules originate in star- forming 
regions and are observed as maser emission in 
two types of sources: class I and class II (e.g., 
Mcntcn 1991). Class II methanol masers are ra- 
diatively pumped and located in the vicinity of 
young stellar objects (YSOs), whereas class I 
methanol masers are believed to trace distant 
parts of the outflows from YSOs and are colli- 
sionally pumped. The most accurate Afi/^ values 
can be estimated from narrow and symmetric line 
profiles with similar shapes. Such profiles are pro- 
vided in interferometric observations of the class I 
masers located in the vicinity of IRAS 16547-4247 
(G343. 12-0.06) which is a luminous YSO with a 
radio jet (Voronkov et al. 2006, hereafter V06). 
This methanol maser emission consists of a clus- 



ter of six spots spread over an area of 30 arcsec 
in extent. One spot (called B in V06) shows ac- 
tivity in 12 CH 3 OH transitions (Fig. 3 of V06). 
Among them, the 9.9 GHz and 104 GHz lines 
have a narrow spike (FWHM < 30 m s _1 ) on 
top of a broader (FWHM w 300 m s" 1 ) sym- 
metric line. The widths of these spikes do not 
exceed the channel spacing which is 29 m s _1 
and 22 m s" 1 at 9.9 GHz and 104 GHz, respec- 
tively. According to V06, these are the narrowest 
spectral features ever found. Their linewidths 
imply that the masers are unsaturated and that 
the turbulent motion in the gas is strongly sup- 
pressed. The 9.9 GHz and 104 GHz masers show 
similar angular sizes, 8 9 9 = 0.10 ± 0.09 arcsec, 
and 9 wi = 0.2 ± 0.1 arcsec, and almost the 
same coordinates, |Aa| = 0.03 ± 0.04 arcsec, and 
\A5\ — 0.09 + 0.04 arcsec. The lines were observed 
in a month interval on June 16, 2005 (9.9 GHz) 
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Fig. 1. — Comparison of the sensitivity coefficients 
calculated in the present paper (Q M ) with those 
from Jll (—Kfj). Shown by error bars are the 2a 
uncertainties. 



and on Aug 18, 2005 (104 GHz). 

As noted in V06, the 9.9 GHz and 104 GHz 
transitions belong to the same J_2 — {J + l)-i E 
transition series with J = 8 and 10, respectively, 
and are expected to show a similar behavior, i.e., 
both maser lines should originate in the same vol- 
ume and, hence, have equal radial velocities. The 
radial velocities of these spikes, V9.9 = —31.554 
km s _1 and Viqa = —31.594 km s _1 , are measured 
with the uncertainty of a few m s _1 , but, unfor- 
tunately, their difference, AV — V9.9 — V104 = 40 
m s _1 , is less certain because of the errors in the 
rest frequencies, £9.9 = 60 m s _1 and £104 = 20 
m s _1 (see Table [TJ. Since at these two fre- 
quencies the sensitivity coefficients are slightly 
different for two calculations [Q^.g.g = —14(1), 
-K^g.g = -11.5(6), and Q M , 104 = -0.45(16), 
—^,104 = —0.18(5)], we use their average val- 
ues Q = (Qf_t — K IM )/2 to estimate Afi/fi from 
Eq.rjU Thus, with Q M ,9. 9 = -12.75(58), Q^ wi = 
-0.32(8) and AV = 40(63) m s" 1 we find Afi/fi 
= 11 ± 17 ppb (lcr, c.l.) or the upper limit 
A/1//1 < 28 ppb. This single point estimate 
is in line with a limit derived from the sample 
mean of ammonia observations in the Milky Way, 



An/n = 26 ± 3 ppb (lcr, c.l.), or \Afi/fi\ < 29 ppb 
(Levshakov et al. 2010a). 

The reliability of the present estimate of Afi/fi 
can be further improved if new laboratory fre- 
quencies of the methanol 9.9 GHz and 104 GHz 
transitions will be determined with a higher accu- 
racy. Since maser sources are in general variable 
in time, additional gains can be obtained from si- 
multaneous observations of several methanol lines 
in a way described, e.g., in Voronkov et al. (2011) 
where up to eight methanol maser transitions were 
observed simultaneously with the Australia Tele- 
scope Compact Array (ATCA). Unfortunately, the 
observed methanol profiles at 24, 25 GHz and 
9.9 GHz have close sensitivity coefficients with 
AQ w 3 which is comparable to the ammonia 
method. This does not allow us to improve the 
aforementioned upper limit on A/x//i. Without 
new laboratory studies and improvements in as- 
tronomical observations, which require substan- 
tial care to determine frequencies with an accu- 
racy better than 10 -8 , any further advances in 
exploring Afj,/ fx from methanol maser spectra will 
be impossible. 
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G343. 12-0.06, Irina Agafonova for her comments 
on the manuscript, and Christian Henkel for many 
helpful discussions. The work has been supported 
by the grant No. 'SFB 676 Teilprojekt C4', the 
RFBR grants No. 09-02-00352 and No. 11-02- 
00943, and by the State Program 'Leading Sci- 
entific Schools of Russian Federation' (grant NSh- 
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